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Abstract 

A basic problem of approximation theory, the approximation of 
functions from the Sobolev space Wp ([0, l] d ) in the norm of L q ([Q, l] d ), 
is considered from the point of view of quantum computation. We 
determine the quantum query complexity of this problem (up to log- 
arithmic factors). It turns out that in certain regions of the domain 
of parameters p, q, r, d quantum computation can reach a speedup of 
roughly squaring the rate of convergence of classical deterministic or 
randomized approximation methods. There are other regions were the 
best possible rates coincide for all three settings. 



1 Introduction 

We are concerned with the study of numerical problems of analysis in the 
quantum model of computation. A series of papers dealt with scalar valued 
problems, that is, with problems, whose solution is a single number. In [H] 
for the first time vector (function) valued problems were considered. The 
first analysis for such a type of problem with matching upper and lower 
bounds was carried out in [H] . 

The present paper is a continuation of [§]. We study one of the basic 
problems of approximation theory, the approximation of functions from the 
Sobolev class W;{[0,l} d ) in the norm of L q ([0, l] d ), a problem which has 
received much attention in the past in the classical settings (see the survey 
[2] and references therein). 

Our results show that for p < q, the quantum model of computation 
can bring a speedup roughly up to a squaring of the rate in the classical 
(deterministic or randomized) setting. On the other hand, for p > q, the 
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optimal rate is the same for all three settings, so in these cases there is no 
speedup of the rate by quantum computation. 

Our method of analyzing the function approximation problem is similar 
to the one developed in [5], namely, we discretize the Sobolev embedding 
problem and show that a sufficiently precise knowledge about the discrete 
building blocks, the embeddings of finite dimensional L p into spaces, 
leads to a full understanding of the infinite dimensional problem. Although 
in a completely different setting, this is close in spirit to Maiorov's discretiza- 
tion technique from approximation theory These finite dimensional em- 
beddings were studied in 0, the results of which will be exploited here. In 
this sense the present paper is related to in a similar way, as a previous 
paper [Sj on quantum integration in Sobolev spaces was related to results 
on summation jHIZj- 

For an introduction and notation concerning the quantum setting of 
information-based complexity we refer to Section 2 of Some general 
results which we will use can be found in Section 3 of that paper. Finally, 
we also refer to jH] for comments on the bibliography. In Section 2 of the 
present paper, which contains the main result, we study approximation of 
the embeddings of Sobolev classes WT([0, l] d ) into L q ([0, l] d ). In Section 3 
we shortly discuss the cost of the algorithm in the bit model and compare 
the results to the classical settings. 

2 Approximation of Sobolev Embeddings 

Let D = [0, l] d be the d-dimensional unit cube, let C(D) denote the space of 
continuous functions on D, endowed with the supremum norm. For 1 < p < 
oo, let L P (D) be the space of real-valued p-integrable functions, equipped 
with the usual norm 
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if p < oo, and 

ll/IUoop) = ess su Pte,D I /(*)!• 



The Sobolev space W p {D) consists of all functions / £ L p (D) such that for 
all a = . . . , ay) G A^q with \a\ := a j — r i the generalized partial 

derivative d a f belongs to L p {D). The norm on Wjj(D) is defined as 
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if p < oo, and 

WfWw^(D) = max||0 Q /|| Loo(D) . 

|a|<r 

We always assume that r/d > 1/p. By the Sobolev embedding theorem 
(see E3]), functions from W p (D) are continuous, and therefore func- 
tion values are well-defined. Let B(W£(D)) be the unit ball of the space 
Wp(D) and J pq : W^(D) — > L q (D) the embedding operator J pq f = f (/ G 
WJ(L>)). 

Now we present the main result of this paper. To emphasize the essential 
parts of the estimates we introduce the following notation. For functions 
a,b : N — > [0, oo), we write a(n) Xi og b(n) if there are constants C\,cq, > 0, 
no G N, ai, ct2 G R such that 

ci(log(n + l)) Ql 6(n) < o(n) < c 2 (log(n + l)) Q2 6(n) 

for all n G N with n > uq. Throughout the paper log means log 2 . Fur- 
thermore, we often use the same symbol c, ci, . . . for possibly different pos- 
itive constants (also when they appear in a sequence of relations). These 
constants are either absolute or may depend only on p, q, r and d - in all 
statements of lemmas, propositions, etc. this is precisely described anyway 
by the order of the quantifiers. 

Theorem 1. Let r, d G N, 1 < p, q < oo and assume r/d > 1/p. Then for 
r/d > 2/p - 2/q 

e«(J pq ,B(W£(D)))^ og n- r / d , 
while for r/d < 2/p — 2/q 

el(J pq ,B(W r p (D))) x log n -ar/Wp-a/ff 

Here en(J pq , B(Wp(D))) is the n-th minimal query error, that is, the 
minimal possible error among all quantum algorithms that use at most n 
query calls to approximate J pq on B(W p (D)), in the norm of L q (D) (see (Hj 
for the definition of en). Theorem ^ is a direct consequence of Propositions 
Q and |21 which are stated and proved below and which also contain the 
logarithmic factors. First we derive the upper bounds. 

Proposition 1. Let r, d G N, 1 < p, q < oo and assume r/d > 1/p. Then 
there exists a constant c > such that for all n G N with n > 2 the following 
hold: For p < q and r/d > 2/p — 2/q, 

eUJ Pq ,B(W;(D))) < cn^/>gn) 2 /P- 2 /*, (1) 
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for p < q and r/d = 2/p — 2/q, 

eUJ pq ,B(W;(D))) < cn- r / d (logn) 4 /f- 4 ^ +1 (loglogn) 2 /P- 2 ^, (2) 
for p < q and r/d < 2/p — 2/q, 

eUJ pq ,B(W;(D))) < cn -2r/d+2/ P -2/ 3) (3) 

and for p > q, 

eUJ Pq ,B(W;(D)))<cn- r / d . (4) 

Proof. We need some preparations. We show that the discretization tech- 
nique developed in [5], properly adapted, works also for the approximation 
problem. For the sake of completeness, we recall also needed details from 
0. For I G N let 



2<m_ 1 



f(2 l (s-s u )) if s G D u 
otherwise. 



i=0 

be the partition of D into 2 dl congruent cubes of disjoint interior. Let su 
denote the point in Du with the smallest Euclidean norm. Introduce the 
following operators En and Rn from !F(D,T&), the set of all real-valued 
functions on D, to J~(D, R), by setting for / G J-(D, R) and s G D 

{E H f){s) = f{s u + 2- l s) 

and 

(Ruf)(s) = 

Let P be any operator from C{D) to L ao (D) of the form 

3=0 

with tj G D and ip» G L 00 (D). Assume furthermore that P is the identity 
on Vr-i(D), that is, 

P/ = / for all f€V r -l(D), (5) 

where V r -i(D) denotes the space of polynomials on D of degree not ex- 
ceeding r — 1. (For example, for d = 1 one can take Lagrange interpolation 
of appropriate degree and for d > 1 its tensor product.) Since r > c£/p, 



4 



we have, by the Sobolev embedding theorem P, ^U], W£(D) C C(D) and 
there is a constant c > such that for each / G WT(-D) 

\\f\\c(D) < c||/||wj(£>)- ( 6 ) 

It follows that 

re— 1 re— 1 

||-P/||i,(D) < ^2\f(tj)\\\<Pj\\L q (D) <^2hj\\L q (D)\\f\\c(D) <c\\f\\ W r(D) 
3=0 j=0 

(7) 

(in what follows the operator P will be fixed, hence Sj=o ll^j IIl 9 (d) can be 
considered as a constant). For / G Wp(D) we denote 

i/p 

\]a\=r 
if p < oo, and 

\f\r,oo,D = max ||5 a /|| Loo(D) . 

|a|=r 

Now we use Theorem 3.1.1 in 2]: there is a constant c > such that for 
all / G W;(£>) 

i nf 11^ ~ ^llwjfD) < c\f\ r ,p,D- (8) 

By ©, (0) and ©, 

< inf r Mf-g)-P(f-9)\\L qi D) 

- c j ni ,„A\f ~ 9\\wr(D) <c\f\ r ,p,D- (9) 
geVr-i(D) 

For I G N set 

2 di -l 2 d! -lK-l 
i=0 i=0 j=0 

Note that 

11^/11^(^=2-^11/11^(0) (feL p (D)). (10) 
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Then we have for u G {p, q} and all / G W£(D), using @ and (fTU|). 

2 di -l 

11/ - ^/||i„(£)) = II E (RljElif — RliPEuf)\\L u (D) 

i=0 

E ll(^(^/- p ^/)IILp) 

8=0 

2 di -l x 7 



2-* e ii%/-^h/iil 

1/m 



(0) 

i=0 



2 a '-l 



< c 2- dl E \Elif\ U r, P 

\ i=0 

< c2 max(l/p-l/n,0)d« | g-dJ |-%/|£ p £, 



2<"-l * 1/p 



and 



j=0 



z- dl E I^k/I^J = (jr* E E J D \d a f(su + 2- l tWdt 



i=0 / \ i=0 |a| 

/ 2 di -l 



2_H E E i^/wp* 



i=0 |a|=r' 

= 2- ri |/|r, P ,D< 2- ri ||/|| W;(D) (11) 

(with the usual modifications for u = oo or p = oo). Consequently, 
ll/-^/lkp) < c2- w + max ( 1 /P- 1 /«'°) d Vk P ,D 

< c 2-ri+max(l/p-l/«,0)«fl||y|| wy(£))i ( 12 ) 

Similarly to [3], we first approximate / by P[*f for some /*, giving the 
desired precision, but using a number of function values much larger than 
n. This Pi*, in turn, will be split into the sum of a single operator P/ , with 
number of function values of the order n, which we compute classically, 
and a hierarchy of operators P[ (I = Iq,...,1* — 1). We will show that 
the approximation of the reduces to the approximation of appropriately 
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scaled embedding operators Jp g l : L^ 1 — > L^ 1 for suitable N\. This enables 
us to apply the results of jS]- Define 

P'f := (Pi - P )f 

2 d -lK-l K-l 

= EE^M+^ffi-E/feiw. ( 13 ) 

i=o j=o i=o 
which can be written as 

p, / = E (E a i*/W*) 

jr=0 V fc=0 / 

with 

k" < «(2 d + 1), (14) 

cijk G R, tj k £ D (j = 0, . . . , k' — 1, k = 0, . . . , k" — 1), and a linearly 

independent system (ipj)jJ^ C L oc (D). The linear independence implies 
that for ii G {p, (?} there are constants ci , C2 > such that for all ay € R 
(j = 0,...,k'-1) 

cilKai^lLg' < II E "iV'ilkp) < ^HK^IIls'- (is) 

Put 

V'Zij = -Rji^j, (16) 

and let 

H = span{Viy : « = 0,...,2 rfJ -l, j = 0, . . . , k' - l} C L 00 (D). 

By the disjointness of the interiors of the Dn and by (|1(J|> we have for ay € R 
(i = 0, . . . , 2 d ' - 1, j = 0, . . . , k 1 - 1) and 1 < u < oo 

2 dl -l k'-X 

II E E a y^yHL(o) 

2 di -l k'-1 2 di -l re'-l 

= E II E a ^'^illL(B) = E E ^i^'IIL(D) 

i=0 j=0 i=0 j=0 

2 d '-l re'-l 

= 2 ~ dl E iiE^-iiLw ( 17 ) 

i=0 j=0 
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Let Ni = K'2 dl . Then 

2 di -l 

E IIKOto 1 !!^ =IIK0II>p (is) 

8=0 

where (aij) stands for (QfyOi=oj=o -1 - Combining l|17|). ((TH|) and (|T3|) . we get 

ci||(ay)ll^ ^ II E E^'^'Hm^) < c 2||(aii)ll L ^- (19) 

i=0 j=0 

Relation ()19|) holds also for u = oo, which can be proved with the usual 
modifications in the reasoning above. Define the operator 7} : II; — > by 

2 di -l re'-l 

J| XI E "ij'V'Jy = (2°) 
i=0 j=0 

It follows from (JTUJ) that for / E II/, 

ll^/H^ <c^\\f\\ Lu{D) , (21) 

and for g G , 

\\Tr l g\\ Lu{D ) <P2NI L ^. (22) 
For Z G N and / G C(L>) set 

ifi/ = RuP'Euf =EE fl J* + r Vi'i- (2 3 ) 

j=0 fc=0 

2 di — 1 

= E P «" (2 4 ) 
i=0 

It is readily verified that 



2«_i 



p i+ i = e ^a^. 

and therefore 



i=0 



2 di_ 1 

-Pz+i — -fi = E Ru(PiEu — PqEu) 
1=0 
2 d '-l 

= E p «= p /- ( 25 ) 

i=0 



S 



From ([T0|). (JI2J) with u = p and lfTT|). we get 

\\ P lf\\L p (D) 



j=0 



2<«-l ^ VP 



2"* £ ||Pi^/-Po^/|| P ip(D) 

' 2<«-l ^ 1 ■'' 

< I 2~ d < Yl (W E l*f - PiEiif\\ Lp (D) + \\Euf - P Euf\\ Lp{D) ) P 

\ i=0 

< cW dl \Elif\ P r>p>D J <c2- rl \\f\\ W r (D) . (26) 

We define operators XJ\ : W£(D) — > by 

17, = T,P/ (27) 

and V; : Lj*« -» L q (D) by 

^=Tf 1 . (28) 

Then clearly 

ViJ^Ui = P(, (29) 
moreover, by ((2*67) an d IIP f° r u = P> 

||^/||^<c2- ri ||/|| w , (I3) (/eWJ(D)) (30) 

and, by (|2*2*|) for u = q 

ll%[U,P) < clblU G?eLfO- (31) 



Now we are ready to derive the upper bounds. It obviously suffices to 
prove them for 

n > max(K, 5). (32) 

Define 

h=ldT 1 log(n/K)\. (33) 
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Then Iq > 0. Furthermore, let 

r = { 1°, * p ~ q (34) 

\ 2l if p<q. 

By & 

^ = Pi + E p /- ( 35 ) 

l=lo 

In the sequel we consider the Pi and P/ as operators from Wp(D) to L q (D). 
By (J33J), K2 d '° < n, hence 

eS(n,,B(^(Z>)),0)=0 (36) 

(this just means that with k 2 dl ° queries we can compute Pi , classically, or, 
more precisely, up to any precision by simulating the classical computation 
on a suitable number of qubits). Let fi,ni € N (I = Iq,...,1* — 1) be 
natural numbers which will be specified later on, and which will be assumed 
to satisfy 

I>-^ 8 <~. (37) 

l=lo 

Put 



J*-l 

n = n + 2k" E v l n l (38) 
i=i 

(if /* = Iq, we do not define the numbers ui,rii and put n = n). From (j!2|) 
above with u = q and Lemma 6(i) of [H, we get 

el{J pq ,B(W; (£>))) 

< sup \\J pq f - I\.f\\ Lq( D) + el{P v *,B{W;{D))) 
feB(W-(D)) 

< c2 -rr+max(l/p-l/5,0)dr + ^(fl, , g(Wj(.D))). ( 39 ) 

The upper bound for the case p > q follows directly from (|39|) and (|36j) . since 
in this case I* = Iq and n = n (this is the trivial case where the optimal rate 
is already attained by a classical algorithm). 

In the rest of the proof we assume p < q. By Lemma 3 of |Hj and Q36JI . 

< em ,B(W;(D))^) + 4_ n (Pi^-Pi ,B(W;(D))) 

= el_ n (P r -P lo ,B(W;(D))). (40) 
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From flsSJ), (J3HJ)- Corollary 3 of 0, and (J37J) we get 



T-l 



< 2^4 K »n i (^^(% r (^)))- (41) 
l=l 

Using Lemma 2 of 6 , (|29|) . and (|31j) we obtain 

^(^,5(1^(2)))) = e^V^UuBiW^D))) 

< cel K „ ni (J^U u B(W;(D))). (42) 



Joining relations (|39|) - (|42|) . we infer 

4(j pq ,B(w;(D))) 



< c2 ~n* + (i/P-iNW + c £ e^(J#ff,,B(W£(I?))). (43) 

In a further reduction one would like to remove the £7; in the last relation. 
This could be done on the basis of Corollary 1 of jS], if the XJ\ were of the 
required form, which is not the case. Instead, we shall approximate the XJ\ 
by appropriate mappings T t : B{W r p (D)) -> (l < I < I*). Note that 

by (§UJ), (El, (El, and (f27|). 

W, J) = E + 2 ~^)' ( 44 ) 

fc=0 

Fix an m* € N with 

2 -m*/2 < + !)-i 2 -rf* (45) 

and 

2 m*/2-l > C) ^ 

where c is the constant from ©. Hence, 

ll/llc^^S™*/ 2 - 1 for feB(W;(D)). (47) 

Define /? : R -» Z[0, 2 m * ) for z G R by 

if z < -2 m */ 2 - 1 

P(z) = { [2 m */ 2 (z + 2 m */ 2 - 1 )J if - 2 m */2-! < z < 2 m */2-i ( 48 ) 
2 m * - 1 if z > 2 m */2-i 
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and 7 : Z[0, 2 m *) —> R for y G Z[0, 2 m * ) as 

7 ( y ) = 2 " m */2 y _ 2 ™*/2-i_ (49) 

Then we have for -2 m */2-i < z < 2 m */ 2 - 1 , 

7(P(z))<z< 7 (P(z)) + 2' m '/ 2 . (50) 

Define m ■ Z[0, JV,) -> D (A; = 0, . . . , re" — 1) by 

»7zfc(*, J) = % + 2-'^ (0 < i < 2 dl - 1, < j < re' - 1) 

(here Z[0,JV z ) stands for {0,1,..., TVj - 1} and we identify Z[0, Ni) with 
Z[0,2 dl ) x Z[0,k')). Next let gi : Z[0,JVj) x Z[0,2 m *) K " -> R be given by 

Qi{{i,j),Vo, • • • = X] a ikl{Vk)- 

k=0 

Finally, we define Ti : B{W r p {D)) -> by setting 

WXi, j) = «((*,.?), (/3 o / o j))^ 1 ). 

for / e #(W£(.D)). Note that T/ is of the form (4) of needed to apply 
Corollary 1 of that paper, which we will do later on. We have 

k"-1 
k=0 

hence, by pi|). ((171) . J5UJ), and (gS)) 

IW)(i,i)-r,(/)(i,j)| 

< E + 2~%) - j((3(f(su + 2-% k )))\ 

k=0 

k"-1 

< 2' m *' 2 E |ojfc| < c2" m */ 2 < c(l* + l)- 1 2- rf * J (51) 

fc=0 

and therefore, for all / € B(W£(D)) and it G {p, <?}, 

Il^/-r z (/)|| z ^ 

(2 rfZ — 1 re'— 1 
(/s / r i2-* £ El(^/)^^)-rK/)(^i)l 
i=0 j=0 

< c(/* + l)- 1 2- w *. (52) 
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Moreover, by (|3U|) and (|52(1 with u = p, 

\\Tl(f)\\ L »i < \\Uif\\ L N l + \\T l (f)-U l f\\ LNl <c2' rl . 
Consequently, 

Fi(B(Wp(D))) C c2~ rl B(Lp l ). (53) 
From (|52|) with u = q and Lemma 6(i) of [I] 

< c^ + lj-^-^+e^^T^^p))). (54) 
Corollary 1 of relation (|53j) above and Lemma 6 (hi) of [3| give 
4^(^*0^))) < el l (J$,c2-rl B (L»')) 

From (|54|1. (|55|) and (|43[). we conclude 
e3(J M ,B(W; (£»))) 



c2- ri e^(J^,S(L^))- (55) 



i i 



+ c^2-^(J^,^)). (56) 



Thus we reached the desired reduction and can now exploit the results for 
the finite dimensional case: By Proposition 2 of jS] and ()56|) . 

el(J pq ,B(W;(D))) 
< c 2~ rl *+^ dl * 

+cY J ^ rl ni ( ^~ q) Nf"^log(n l /y^ + 2))p-l. (57) 
i=i 

Recall that we consider the case p < q. First we assume 

r 2 2 

- > . (58) 

a p q 



Fix any 5 > such that 



r >[---) (d + 8) (59) 
p q 
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and put for I = Iq 



ni 



^dl -5(l-l ) 



(60) 



vi = [8(21n(/-/ + l) + ln8)l . 
It follows from H61|) that 

E^ i/8 4&- /o + ir2< i' 

l=lo l=lo 
so assumption (|3*Tj) is satisfied. By (|3*5|) . (|6*Tj) . |23|)> and (|3"4l . 



(61) 



(62) 



z*-i 



n 



< n + 2k" ^ [8(2 ln(/ - l + 1) + In 8)] 2 dl °- s( - 1 - 1 ^ 
l=lo 

c2 dl ° < cn. 



It follows from (J57|), (fS^-flMj) . and © that 

e ^(j M ,s(w;(D))) 
i*-i 



(63) 



< c2-^/2 + c ^ 2 - 



W-(^-|)d« +(^-f)5(i-/o)+(f-f)^ 



l=lo 



= c2~ rl ° +c2- rlo {l + l)v— q J^2 

;=z 

< c2~ r/ °(/ + 1)p - « < cn~5(l gn)p~9 



J : (-H-(|-f)(<H^))(J-to) 



(64) 



Then Q follows by an obvious scaling from (|63|) and (|64|) . Next we assume 



r 2 2 
d p q 



Here we take any 5 > with 



r < ( - --) Ul-*) 



and put for I = Iq, . . . ,1* — 1 



?dl -6(l* -/) 



(65) 



(66) 
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vi = [8(21n(r -I) + In 8)1 . 
From (|H7j) we see that (pT7|) is satisfied, again: 

1=Iq l=lo 
By ((SHI), &, ©, JHEJ), and 

r-i 



7? 



< n + 2«"^ [8(21n(i*-0 + In8)l 

i=i 

< c2 dl ° < cn. 

Using JS7|), (EU)-®, and ©, we get 

e 9(j M ,B(w;p))) 

< C 2 ( "^ + I"l )d ' 

- r j_(|-|)«Bo+(2-fW-0+(i-2)dJ 



? dZ -5(«*-0 



(i*-i + iy 



Z=Z 

< c2 d p i' < cn d p ?, 

which, together with (|69jl. gives ©• Finally we consider the case 

r _ 2 2 
d p q 

Here we put for I = Iq, . . . , I* — 1 

2^(Zo + ir 1 (ln(Z + 2)r 1 " 



ni 



and 



vi = [8(ln(Z + 2)+ln4)l . 
This way, relation (|H7|) is valid: 

r-i , r-i 
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From (|3SJ, (J7UJ), (JZU, (E3J), and we get 

r-i 

n < n + 2k" ^ [8(ln(Z + 2) + In 4)] ^ (/ + l)" 1 (ln(J + 2))" 1 

< c2*><cn. (73) 

By the help of (JS7J), ip ^ -IpE ]) . and J7DJ) we conclude that 

el(J pq ,B(W;(D))) 
l*-i 

< c2 - w */2 + c J- 2 - rl - {£ p-^ dlo+ ^-^ dl (l + 1)1-1 (log(Z + 2))t" 



^ ,V 2 -H-(|-|)*+(|-f )* (Zo + l)|-!(l og(Zo + 2))|-i 

r-i 

2 2 ^ — - , 



c2- rl ° + c{l + l)F"(log(fo + 2))"i 2 

< c2- w o(/o + l)t-^ +1 (log(; + 2))l " 

< cn d (log n)p 9 (log log n) 



2 2 
<3 



2 _ 2 
)P 1 , 



This shows © and completes the proof. □ 

Proposition 2. Lei r, d 6 N, 1 < p, g < oo, and suppose r/d > 1/p. Then 
there exists a constant c > such that for all n € N with n > 4 i/ie following 
holds: First assume r/d > 2/p — 2/q. If 2 < q < oo, then 

el{J pq ,B{W;{D)))>cn-rl d , (74) 

if q = 2, then 

eUJ P ,2,B(W;(D))) > cn-^(loglogn)- 3 / 2 (logloglogn)- 1 , (75) 
and if 1 < q < 2, then 

eUJ Pg ,B(W;(D))) > cn-^Oogn)- 2 /^ 1 . (76) 
Now assume r/d < 2/p — 2/q. Then 

el{J M ,B{W;(D))) > cn -2r/d+2/ P -2/ q _ (7?) 

Proof. Let ^ be a C°° function on R d with 

supp^ C (0,l) d , o\ := / ^(t)dt>0, 

Jd 



16 



and denote ||V'|lw r (D) = °"2- Let n € N, G No, and iV = 2 dfc . Let and 
be as defined in the beginning of the proof of Proposition^ Set 

il>i(t) = Rkrt (i = 0,...,N-l). 

We have 



il>i(t)dt = criN- 1 (78) 



and 

Ui\\wr { D) < 2^- d M k U\\ W r (D) = a a 2<r-W k . 

Consequently, taking into account the disjointness of the supports of the ipi, 
forallaiGR (i = 0, . . . , N - 1), 

N-l /N-l \ 1 / p 

II E ^ ™ m = E wi^r™ < ^IK^IU ( 79 ) 



1=0 pV ; \ 1=0 y 



P 



(which holds also for p = oo). Fix any m* 6 N with 

m*/2 - 1 > dk/p. (80) 

Let (3 : R -> Z[0,2 m *) and 7 : Z[0,2 m *) -> R be defined as in (gSJ) and (jMJ). 

For / € B(Lp) we have 

Hence, by ([50]). for < i < N, 

7 WW)) < /(*) < 7W(0)) + 2- m * /2 - (81) 

Define 

JV-l 

r:B(L^)-^Wj(£)) by r(/) = £ 7 ° /? o f(i) fr. 

i=0 

By (jTHJl and (JUJ), for / G B{L$), 

\\T{f)\\ W r {D) < a 2 2 rk hal3of\\ L „ 

< a 2 2 rk (\\f\\ L » + ||/- 70/30 

< a 2 2 rfc f 1 + 2- m */ 2N ) . (82) 
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Define $ : L q {D) -» L*j by 

(*/)(*)= iV / f(t)dt. 

JD kl 

It follows from (J7BJ) that 



= o-i e*, (83) 

-TV 

IC LiCllULCO 111c t-tn unit VCl^tUl 111 A 

ity, 

ZV-1 



where ej denotes the i-th unit vector in . Moreover, by Holder's inequal- 



i=0 ^ D fci 
ZV-l „ 

< ^E / i/ww^r 1 



L q (D)- 



i=0 

Thus, since is linear, 

||$||Lip = \M < 1. (84) 

Furthermore, by (JB3J), 

N-l 

t>°J pq oT(f) = ^ 7 o/3o/(i)^ 
i=0 
ZV-1 



i=0 

= ^(7^°/). (85) 

Define 77 : D -> Z[0, iV) by 

7?(s) = min{f | s G 

and 

g:DxZ[0,2 m VR by g(s, z) = l{z)% {s) {s). 

Then 

N-l 

r(/)(«) = E 7 °' 9o ^)^( s ) 

= 7°^/(l(«))^( s )(s) 
= q(s,P ° f °v(s))- 
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So T is of the needed form (see relation (4) of [5], with k = 1) and, by 1)82(1 . 
maps 

B(L%) into a 2 2 rk (l + 2- m *' 2 ^B{W;{D)). 
By (|84|). Lemma 2 and Corollary 1 of [0], and Lemma 6 (hi) of [3], 

4(<Doj M or,s(Lf)) 

< 4(j M or,B(if)) 

< e£ (j pq ,a 2 2 rk (l + 2- m */2) ^(D))) 

= a 2 2 rk (l + 2- m */ 2 )eUJ Pq ,B(W;(D))). (86) 
Using (|81|) again, we infer 

sup || J^f - J*( 7 o p o /)|U = ||/ - 7 o /? o /|U < 2" m */ 2 , 

and hence, by Lemma 6(i) and (ii) of 4_, and (JBSJ), 

el n {J^B{L»)) < e^(J^oW,^)) + 2— */2 

= a 1 - 1 eL(<I>oJ M or,e(L^)) + 2— V2 

< arV 2 2 f ' fc (l + 2- m */2^) e q ( Jp9) (£»))) + 2" m */ 2 

< d\r/ d ea(J M) B(W£(£>))) + 2^ m */ 2 , 
where 70/3 stands for 

(7o/3,..., 7 o/3):R Jv ^R iV . 

Since m* can be made arbitrarily large, we get 

e* n (J»,B(L»)) < cN r / d eUJ Pq ,B(W;(D))). (87) 

For the case r/d > 2/p — 2/q, we choose k = \d~ l (\og{n/ co) + 1)], where Co 
is the constant from Proposition 6 of [5], which can be assumed to satisfy 
< Co < 1. It follows that 

c 2- d N = c 2- d 2 dk <2n< c 2 dk = c N. (88) 

Now the lower bounds (|71|). (|7H|). and JZHJ) follow from (jEZj), (|HH|, and Propo- 
sition 6 of 0| . In the case r/d < 2/p — 2/q, which implies p < q, we set 
A; = |~<i -1 (log(n 2 /co) + 1)] . Consequently, 

c 2~ d iV = co2~ d 2* < 2n 2 < c Q 2 dk = c N. (89) 

Relation (|77j) results from (|87j) . (|89|) . and Proposition 4 in jSJ. □ 
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3 Comments 



The algorithm we presented was optimal with respect to the number of 
queries. (Although parts of the algorithm occur only in an implicit way, 
through the use of properties of e« numbers, it is straightforward to trans- 
form all upper bound proofs into algorithmic details.) Let us now consider 
its cost in the bit model of computation. Here we assume that n and N are 
powers of 2. We use the respective remarks about bit cost made in Section 
5 of 0. 

For p > q classical approximation suffices. For p < q the problem 
splits into the classical computation of Pi and the approximation of J^ 1 
for I = l , . . . , I* — 1 = 2/q — 1 = C(logn) using rii queries (see the proof 
Proposition ^ for these numbers). To increase the respective success prob- 
abilities appropriately, we have to repeat these approximations V[ times on 
level I, and we have V[ = (log log n). The total number of queries is n (or 
n = 0(n) if considered before scaling), that of quantum gates is 

\l=k J 
The algorithm needs O(logn) qubits and 

O i "WlNf 1 max(log(n i /0Vi), l)" 1 J = 0(n/ log n) 

measurements. To compute P[ f classically, we need 0(n) function values 
and O(nlogn) classical bit operations. For the approximations on the levels 
a total of 

O ( Y WiNf 1 log N t \ =0(ri log n). 
\l=lo ) 
classical bit operations is required. This does not yet take into account the 
classical computation of the vector analogue of the median. Let us assume 
that we apply the constructive procedure described after Corollary 1 of [S] . 
At level I we have to compute the norm of vf vectors in L^ 1 with at most 

O (nfNf 1 max(log(V7A,), l)" 2 " 

non-zero coordinates. This amounts to 

/ ^ _ 
O logn ^ i/ i 2 nfA i ~ 1 max(log(n i /v / A^),l)" 2 I =0(nloglogn) 

\ l=lo 
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classical bit operations. We see that the overall quantum bit cost differs by 
at most a logarithmic factor from the quantum query cost O(n). 

The concrete form of the output of the algorithm depends on the struc- 
ture of P. If P is, for example, tensor product Lagrange interpolation, then 
the output is a sum of O(logn) piecewise polynomial functions, with 0(n) 
pieces for the classical part P[ f and 

O (nfN- 1 max(log(n z /v^), I)" 2 ) 
pieces not identical to zero on level I, that is, a total of 

O | n+ ^nfiVf 1 max(log(n i /v / iVi),l) _2 ) = 0{n) 

nontrivial pieces, with each point of D being contained in at most O(logn) 
pieces. 

We summarize the results on the approximation of J pq : B(W£([0, l] d )) — > 
L q ([0, l] d ) in a table and compare them with the respective known quanti- 
ties in the classical deterministic and randomized settings (see jS] and the 
bibliography therein). Recall that we always assume r/d > 1/p. The respec- 
tive entries of the table give the minimal errors, constants and logarithmic 
factors are suppressed. 



Jpq 


deterministic 


random 


quantum 


l < p < q < oo, 


n -r/d+l/p-l/q 


n -r/d+l/p-l/q 


n -r/d 


r/d > 2/p - 2/q 








l < p < q < cxd, 


n -r/d+l/p-l/q 


n -r/d+l/p-l/q 


n -2r/d+2/p-2/q 


r/d < 2/p - 2/9 








1 < 9 < P < 00 


n -r/d 


n -r/d 


n -r/d 



We observe a possible improvement of n~ l (for p = 1, q = 00) over the 
classical deterministic and randomized case (which is essentially a squaring 
of the classical rate for r/d close to 1). This is the maximal speedup over the 
randomized case observed so far in natural numerical problems (the same 
speedup was first found in [Sj for integration of functions from W[(D)). We 
also see that there are regions of the parameter domain where the speedup 
is smaller, and others, where there is no speedup at all. 
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